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SUMMARY 

With the use of the wind-tunnel measurements of Bagnold and Zingg, a model is 
developed for estimating the parameters that describe the flux of sand on Mars. Applica- 
tion of this model to the sensor -breakage problem indicates that the expected lifetime on 
Mars of the wind sensors of the Viking 75 Meteorology Instrument System is about 40 Earth 
years for the worst case considered; this expected lifetime is adequate for both the primary 
Viking 75 mission and for a proposed extended mission of one Martian annual cycle. 

INTRODUCTION 

The Viking 75 Meteorology Instrumentation System contains three "hot -film" 
sensors; each sensor is about 0.5 mm in diameter and 1 cm in length and is mounted as 
shown in figure 1. The sensors are made of borosilicate glass coated with a platinum 
film 0.5 pm thick and overcoated with 5 pm of aluminum oxide. Two of these sensors 
are used to measure two components of wind velocity; the third serves as a reference 
temperature sensor for the wind measurement. These sensors are deployed from the 
Viking lander on a nonadjustable boom. Variation of lander tilt within its nominal range 
results in a sensor -height range above the Martian surface of 1.2 to 1.6 m. While de- 
ployed, the sensors are exposed to impact and possible breakage by windblown Martian 
sand. Environmental criteria of the Mars Engineering Model have been developed for the 
exposure of the Viking lander to windblown sand and dust. However, these criteria were 
developed primarily for use in designing the protection of lander components from surface 
erosion, and may not be appropriate for evaluation of the likelihood of sensor breakage. 

The purpose of this paper is to develop and to apply criteria suitable for the evalua- 
tion of the expected lifetime of the hot -film sensors in the Viking Meteorology Instrument 
System while the sensors are exposed to saltating (blowing) sand after landing on the planet 
Mars. It is the intent of this report to develop realistic rather than ultraconservative val- 
ues. However, the uncertainties involved both in the understanding of the saltation pro- 
cess and in the estimation of Mars parameters require the use of conservative assumptions 
and values in several cases. Therefore, the estimated lifetime is a conservative value. 



although a much less conservative one than the erosion criteria for the Mars engineering 
model. The order of the development is: first, criteria for sensor breakage by sand- 
grain impact are developed; second, a saltation model is developed to permit calculation 
of values of parameters needed in the failure criteria; and third, quantitative values for a 
range of expected Martian conditions are applied to determine the expected sensor life. 

SYMBOLS 
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nondimensional constant in flux equation 
sphere drag coefficient 
empirical constant in equation (6) 
standard diameter in flux equation (250 /jm) 
diameter of sand particle 
kinetic energy (KE on computer printout) 
slip factor 

acceleration due to gravity 
scale height = 
mean height of saltation 

height where slope of logarithmic wind profile changes (Bagnold's "focus") 

von Karman's constant (taken as 0.4) 

mass flux per unit area 

mass flux per unit length 

Reynolds number 

wind speed 
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threshold wind speed at height k' 

V* friction velocity ^Tq/p 

V*' equivalent friction velocity during saltation 

^ threshold friction velocity for saltation 

Wq initial vertical velocity or specific impulse 

z height above planetary surface (Z on computer printout) 

Zq roughness length in equations (9) to (11) 

^max apex height of trajectory 

X molecular mean free path 

p coefficient of viscosity 

p atmospheric density 

CT particle density 

Tq surface shear stress 

Subscripts: 

M Mars 

E Earth 

FAILURE CRITERIA FOR WIND SENSORS 

Tests conducted for the Viking 75 project by Wendel J. Maegley of Martin Marietta 
Aerospace^ have established at least partial criteria for breakage of the Viking wind sen- 
sors by impacting particles. Because sufficient flight -type sensors were not available 

^W. J. Maegley: Effect of Wind Borne Particles on the Survival of the VMIS Wind 
Sensors. Martin Marietta Aerospace Memorandum 8943-73-14, Jan. 19, 1973. (Not 
generally available.) 
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only one actual flight -type sensor was tested. A few early prototypes of the flight sensor 
were used; the largest number of test subjects were uncoated borosilicate glass rods of 
the type used in the flight sensor. The number of tests with each test subject is shown in 
table I. 

The test procedure involved blowing variously sized screened basalt particles 
through a tube onto the test specimens. Films of the test were used to determine particle 
velocities and number of impacts. For most of the tests the particle -flux rates were a 
few orders of magnitude higher than the particle -flux rates predicted herein for saltation 
on Mars. However, for the single flight -type sensor the particle flux was many orders of 
magnitude higher. 

The results of Maegley's tests are shown in figure 2 and in table I. The diagonal 
line in figure 2 between breakage and no-breakage corresponds to a particle kinetic energy 
of 15 juj. Particles with less than 15 pj kinetic energy did not break the test subjects 
even after very long exposures (>100 000 impacts). Thus, it appears that the particle 
kinetic energy constitutes a suitable criterion for the combination of particle mass and 
velocity required to break the sensor. The value of 15 /ij is adopted as the minimum 
kinetic energy for breakage. 

The number of impacts of particles with a kinetic energy of 15 pJ or higher required 
to break the sensor is a much more uncertain figure. As shown in table I, from 25 to 
24 000 impacts at or above the given energy level were required to break the various test 
subjects. Compared to the uncoated rods, the coated prototype sensors required a much 
larger number of impacts before fracture, although the kinetic energy per particle was, 
at most, only slightly larger. In both cases the failures appeared to be progressive, with 
a number of microscopic pits or chips appearing before the sensor actually failed. The 
data in table I show that, even though there was a very wide range in the number of impacts 
required for breakage, the number required for each of the test conditions where more 
than one test subject was used had a surprisingly small range of number of impacts. This 
finding also suggests a progressive failure mechanism with failure occurring only after a 
more or less fixed number of impacts. 

The flight -type sensor had a much thinner but much smoother aluminum oxide over- 
coating than the prototype sensors, so the sensors might be either stronger or weaker. 

In the case of the flight -type sensor, which was subjected to a much higher rate of flux, 
extensive damage was not noted before breakage; the actual break occurred when two 
impacts came in rapid succession. Furthermore, motion pictures showed that the break 
occurred at a point on the sensor well out from the point of impact. These circumstances 
suggest the possibility of a resonance effect resulting from the two impacts in rapid suc- 
cession. As can be seen in the "Results and Discussion" section, the probability of such 
rapid impacts on Mars would be extremely small. Since the failure of the flight -type 


4 



sensor probably occurred at an anomalously low number of impacts (1100), this number 
offers a conservative, perhaps a very conservative, breakage criterion. Thus, for the 
purpose of the present study, a value of 1000 impacts is adopted for the number of 15 pj 
impacts required to fracture a sensor. 

The data for a mean particle diameter of 1200 pm show that at the higher kinetic 
energy of 43 pj, the number of impacts required to break. the test subject was consistently 
lower than the number required where the kinetic energy was in the range of about 15 pJ 
to 25 pJ. Since the number required at the 43 -pJ level was about an order of magnitude 
less than for similar test subjects at the lower energy values, a criterion of 100 rather 
than 1000 impacts is adopted for a kinetic energy of 45 pJ. 

For particles with a kinetic energy much above that in the tests, the required number 
of impacts would be expected to be much smaller. For some still larger value of kinetic 
energy the sensors would break on the first impact. 

In order to account for the effect of possible impacts with kinetic energy greater 
than the test values, the two points (15 pJ, 1000 impacts and 45 pJ, 100 impacts) are con- 
nected with an exponential curve, giving the number of impacts to breakage as 

pk(MJ)-15 

1000 X 10 L 30 _ 


An exponential curve is used because a straight line would predict breakage at zero im- 
pacts with a kinetic energy of only 48 pJ, a prediction which is clearly unrealistic; the 
other simple curve which might be used (a power law) would be less conservative than the 
exponential curve. (In the calculations presented in "Results and Discussion" for Martian 
conditions, the critical kinetic -energy values are very near the experimental data points, 
so that the precise shape of the curve is not crucial in this particular application.) 


In actual exposure each particle size has a different value of kinetic energy; a crite- 
rion which varies with kinetic energy is not convement to use. Instead, an effective num- 
ber of impacts is defined equal to the actual number multiplied by 


Ek(pJ)-15 


10 


L 30 J 


and a single criterion of 1000 effective impacts is used. 

Very high particle velocities are required for sensor breakage by small particles 
while larger particles can break the sensors at much smaller velocities. (See fig. 2.) 
This finding suggests that the larger particles may be of great importance in the sensor- 
breakage problem. However, in the erosion criteria for the Mars engineering model, the 
larger particles are of minimal importance compared to particles in the 100- to 200-jUm 
range. 
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SALTATION MODEL 


Particles larger than about 50 jum fall rapidly enough in the atmospheres of both 
Earth and Mars so that they cannot be considered fully suspended particles. In high 
winds, particles whose diameter is below about 1000 /um move by a series of jumps in a 
process called saltation. The saltation process has been studied by investigators with a 
wide variety of interests, and semiempirical laws have been formulated to describe some 
aspects of the saltation process. Such laws have dealt most notably with threshold veloc- 
ities required to initiate or to maintain the saltation process and with the total (height- 
integrated) flux of material across a unit width. However, little information is available 
regarding the variations of particle concentration or flux with height, and no information 
has been found regarding the variation of particle velocities with particle size and height 
above the surface. Nevertheless, the latter quantities must be specified in order to deter- 
mine the effect of the saltating particles on exposed components of the landed spacecraft. 

In this section methods of specifying the necessary quantities are described. For 
parameters other than the variations with height of particle size, velocity, and flux, the 
methods and assumptions used in the present study are the same as those used in formu- 
lating the blowing-sand environmental criteria in the Viking Mars Engineering Model, 
M-75-125-3, Revision A. 


Integrated Mass Flux 

A number of measurements have been made of the height -integrated mass flux q 

where 



and Q is the mass flux at height z. Bagnold (ref. 1) suggests the form 


q = 




( 1 ) 


( 2 ) 


where C is an empirical coefficient which depends on the degree of uniformity of grain 
diameter in the sand and D is the predominant diameter of the sand. This equation has 
been found to fit not only Bagnold' s experimental data but also that of subsequent inves- 
tigators quite well. Zingg (ref. 2) found slightly lower values of C and a three -fourths 
rather than a one -half power dependence on diameter. The latter difference is of little 
importance for sand whose predominant diameter is near the standard diameter D of 
250 pm. 
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For the diameter distribution of dune sand given in the Viking Mars Engineering 
Model, with a predominant diameter of about 300 pm, a C of 1.8 corresponding to natu- 
rally graded dune sand, a g of 3.7 m/s^, a p of 1.2 x lO"^ kg/m^, and a corre- 

sponding to a wind speed of 40 m/s at a height of 1 m, equation (1) yields a value of 
q = 0.1 kg/m-s, the value adopted for the Mars Engineering Model. Values used in the 
present study vary with the assumed Mars conditions, but are close to the above value of q. 

Mass -Flux Variation With Height 

Relatively few measurements have been made of the variation with height of mass 
flux (through a unit area), but Bagnold (ref. 1) presents the variation with height for dune 
sand. Zingg (ref. 2) gives the variation with height for several narrow ranges of sand 
diameter. In both cases the variations are roughly exponential with height. 

The variation of mass flux with height is assumed to be described by 



As required by equation (1) the integral of Q is equal to q. The parameter H deter- 
mines the logarithmic rate of decrease of Q with height, and by analogy with atmospheric 
pressure and density scale heights it is called the mass -flux scale height, or simply scale 
height. Setting the partial derivative of Q with respect to H in equation (2) equal to 
zero shows that the maximum flux at a given height z occurs when H = z. 

For the Mars Engineering Model, the mass flux (0.2 exp -2z kg/m^-s) was specified 
such that it produced the largest mass flux in the height region between 0.5 m and 2 m 
where the equipment is mounted on the lander body. The value in the Mars Engineering 
Model is the same as the value given by equation (3) with H = z = 0.5 m. 

Values of the mass flux at z = 0.5 m with q = 0.1 kg/m-s resulting from given 
values of H are shown in figure 3. Specifying a value of H larger than 0.5 results in 
a smaller value of Q; a smaller value of H yields a much smaller value of Q. 

At heights above 0.5 m, setting the value of H equal to 0.5 m gives mass fluxes 
somewhat lower than the maximum mass fluxes. However, as shown in figures 4 and 5, 
at z = 1.2 m and z = 1.6 m (the lowest and highest expected heights of the wind sen- 
sors) the indicated values of flux are still much larger than would be given by smaller 
values of H such as are derived in later sections. Hence, the indicated values of flux 
may be overconservative. 

The scale height H can be related to the "average height of saltation" h (called 
yg^ in Zingg' s notation) as measured by Zingg (ref. 2). Since by definition half of the flux 
takes place above and half below h, then 
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( 4 ) 


so that 



H 



q 


H = -i-h = 1.44h (5) 

In 2 

Zingg (ref. 2) made wind-tunnel measurements of h for various particle sizes and 
wind velocities under Earth surface conditions of gravity and atmospheric density. The 
results are shown in figure 6. Zingg found that his data could be fitted by an equation of 
the form 

h= Czd^/^p^/^V^ (6) 

where the value of is 2350 

If a way can be found to relate the value of h on Mars to the value of h on Earth, 
equation (6) can be used to find the value of h, equation (5) can be used to find H, and 
equations (2) and (3) can be used to find the value of Q as a function of height on Mars. 
This value of Q can then be used with particle velocities derived by numerical trajectory 
simulation and a particle -size -distribution model to derive concentration, kinetic -energy 
flux, particle flux, rate of sensor impacts, and, finally, an estimate of the expected life of 
the sensors. 


Relation of Saltation Scale Heights for Earth and Mars 

The initial vertical particle velocity Wq is assumed herein to be caused by some 
combination of lift, drag, and/or kinetic energy from a previous "jump." All these quan- 
tities may be expected to be functions of the stagnation pressure i- pV^. Since V at 

a 

any particular height is approximately proportional to V,^', it is reasonable to assume 
that 


Vo,M^ 

where a is assumed equal for Earth and Mars. 

For the same Wq, the scale height and the average saltation height on Mars and 
Earth are assumed proportional to the respective apex heights. Thus, 
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(8) 


Hm 

^M _ j 

_^max(^o)]M 

He ' 

1 

^max(^o)Jg 


Particle -Size Distribution 

Since the value of h and, therefore, of H for a given Wg is a function of particle 
size, the variation of mass flux with height is different for each particle size, and the mass 
flux at a given height must be found by combining the mass fluxes of particles of various 
diameters. Thus it is necessary to specify the size distribution of particles involved in 
the saltation processes. It is here assumed that all particles which contribute significantly 
to the composition of surface material are involved, so that the size distribution of parti- 
cles in the surface material is used. 

The distribution of particles by particle size was taken from the dune -sand model 
distribution of the Mars Engineering Model. The proportion for each lOO-jum-diameter 
class IS shown in table II. In forming this table intermediate particles were assigned to 
the even 100-p.m-diameter class whose mass was nearest the mass of the intermediate 
particle. Particles with diameters smaller than 100 fim or larger than 1000 /im are not 
considered separately because of the rapid decrease of number density with particle diam- 
eter in the dune -sand distribution. However, all smaller particles are included in the 
100-/Lim class and all larger particles in the 1000-/ini class. 

Trajectory Simulations 

In order to determine particle velocities under Martian conditions, a series of numer- 
ical particle -trajectory simulations is performed. In the numerical trajectory simulation 
the sand-particle motion is calculated at constant acceleration over each small time step 
using the wind-velocity profile and the slip-corrected drag coefficient described herein. 
Then, using the new particle position and velocity at the end of the time step, a new wind 
speed at the new height and then a new drag force are calculated; the process is repeated 
for the next time step. The process is straightforward once all the necessary parameters 
are defined. 

In these simulations the particle is started out at zero height with some assigned 
initial vertical velocity Wg. This initial vertical velocity can be thought of as represent- 
ing the specific impulse of the particle imparted by some combination of impact dynamics 
and aerodynamic lift due to wind shear and/or particle rotation. In an actual particle tra- 
jectory, the lift might be distributed over some finite vertical height. However, because 
of the concentration of wind shear in the very lowest levels, the effects of lift can be 
approximated by an instantaneously acquired velocity. Aerodynamic lift is not otherwise 
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included in the calculation; particle -particle collisions are not considered as such colli- 
sions would occur too rarely to play a significant role. 

All particles are assumed to be spherical in shape and to have a density of 
2650 kg/m^. Sphere drag coefficients as a function of Reynolds number are calculated 
by a curve fit taken from Morsi and Alexander (ref. 3). This curve fit utilizes a series 
of quadratics in (1/R). The resulting curve for Cq as a function of R is shown in 
figure 7. 


The particle diameters involved are not sufficiently large compared to the Martian 
atmospheric mean free path (about 10 fim) to insure continuum flow. Therefore, a slip- 
flow correction taken from Davies (ref. 4) has been applied to the drag coefficient. The 
slip factor F = 1 + ^[2.514 + 0.8 exp(-0.55d/X)] was designed as a multiplicative cor- 


rection factor applied to the Stokes’-law velocity. Because the present range of Reynolds 
numbers extends beyond the region of applicability of Stokes' law, the correction is applied 
instead to the drag coefficient Cj) described therein. Thus, 


Cp,(continuum) 
CD(slip) = 


The slip correction ranges from about 3 percent for 1000-jum particles to about 
30 percent for 100-jum particles. 

Wind-velocity profile . - When the atmospheric stability is near neutral, the vertical 
profile of the wind speed in the boundary layer can be approximated by a logarithmic pro- 
file of the form 


V(z) = ^ In |- (9) 

ko ^o 

Once saltation is fully developed the atmosphere is well mixed and the atmospheric sta- 
bility is near neutral. However, the wind profile is modified by interaction with the sal- 
tating particles. Bagnold (ref. 1) and Zingg (ref. 2) found that the modified profile could 
be approximated by a combination of two logarithmic profiles; the wind speed up to 
height z = k' is given by 

V(z) = ^^ln|- (10) 

^o o 

and the profile from z = k' to the top of the boundary layer (which is always above the 
top of the lander) is given by 
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( 11 ) 



where 


Vt = V(k’) = ^ln^ (12) 

ko Zq 

(The velocity V(z) is assumed to be zero for z < Zq.) For Earth conditions, Bagnold 
found k* to be about 3 mm and Zingg found k' to vary with particle diameter and to be 
equal to lOd (for an assumed mean diameter of Martian sand of 300 jim, Zingg's relation- 
ship also gives k' = 3 mm). In view of the lower gravity (3.7 m/s2) and possible greater 
saltation height of Mars, calculations have also been made using k' = 1 cm. Values used 
for the roughness parameter Zq are 300 fim corresponding to observed values for Earth's 
deserts and 30 jum corresponding to smooth sand. 

Atmospheric parameters .- The values required for the calculations are the density 
and coefficient of viscosity. Atmospheric density enters into the calculations in several 
ways: it appears directly in the horizontal- and vertical-acceleration terms; it enters 
the drag coefficient through the Reynolds number; and it enters the slip -flow correction 
through the mean free path. Viscosity enters only through the Reynolds number. 

Representative Particle Trajectory 

If the probability -density function of Wq for each particle diameter were known, 
particle velocities (and, indeed, all the variables of interest) could be determined by 
straightforward Monte Carlo simulation. Lacking knowledge of this function a single rep- 
resentative trajectory for a class of particles may be established (i.e., the class of par- 
ticles of a particular size at a particular height) . 

In the present model most of the particles at a given height are particles which have 
apex heights at or slightly above the given height. Particles which have apex heights below 
the given height cannot, of course, collide with a sensor at that height. With an exponential 
decrease of mass flux with height, the mass flux decreases rapidly with height. With in- 
creasing particle velocity with height, the number of particles at a given height decreases 
even more rapidly than does the mass flux. In "Results and Discussion” computed scale 
heights indicate that the decrease is very rapid on Mars. In addition, since the vertical 
velocity of the particle is very low near the apex height, the probability that a particle 
would be at some height below its apex height, and thus could strike a sensor at the lower 
height, decreases rapidly with decreasing height. For example, the calculated probability 
density -per -unit height for a 300 -/Lcm particle with a 2-m apex is shown in figure 8. The 
decrease with decreasing height in the region just below the apex height is very rapid. 

Because of the combination of the rapid decrease with height of the number of parti- 
cles which reach levels above the sensor height, and because of the very rapid decrease in 
the probability that one of these particles might strike the sensors, particles with apex 
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heights equal to the sensor height are considered to be representative of particles which 
strike the sensors. The velocity of such particles is used in subsequent flux calculations. 

Values of particle velocity for various particle sizes and apex heights are deter- 
mined by the mathematical simulation described in the section on "Trajectory Simula- 
tions." The Wq values are selected to produce the desired apex height. The resulting 
particle velocities are shown in table in, and the corresponding particle kinetic energy is 
shown in table IV. The kinetic energy of the smaller particles is less than the 15 pj 
which Maegley found was required for sensor breakage. 

The kinetic energy of the small particles is below that required for sensor breakage. 
This fact is important in the use of the present model for two reasons. First, the simu- 
lations by Zingg in a wind tunnel with a 91 -cm (3 -ft) square test section may not have 
adequately simulated the effects of larger-scale atmospheric turbulence on the smaller 
particles which may be in partial suspension. Sharp (ref. 5) presents data on the "height 
at which each grain -size fraction attains its maximum weight percentage" for measure- 
ments made in the atmosphere of Earth. These data show values of this "height" for 
particle -size ranges of <62 pm and of 62 pm to 125 pm to be larger than the values for 
larger particle-size ranges. For diameters larger than 125 pm the values increase 
monotonically with increasing diameter in general agreement with Zingg' s wind-tunnel 
data. Since fall velocities are about the same on Earth and Mars for particles between 
10 and 1000 pm, similar behavior might be expected on Mars. Second, on Mars there may 
be a much larger percentage of very small particles than the percentage in the dune -sand 
model distribution of particle sizes, since the major small -particle sinks of Earth (e.g., 
oceans) are not found on Mars. If either of these conditions should exist, and if these 
small particles were to possess sufficient kinetic energy to break the sensors, then the 
present model could be unconservative for use in sensor -breakage studies. Since the 
small particles do not possess sufficient kinetic energy, the model is not unconservative 
for the present use. 

Because of these two reservations concerning its applicability to small particles, 
it is not recommended that this new model be used to reevaluate the Mars Engineering 
Model criteria for erosion, since erosion was found to be produced largely by particles 
of the smaller sizes. However, the new model is considered to be applicable to the 
problem of sensor breakage where the smaller particles play no part. 

Threshold Velocity 

In the calculations described above it is assumed that the friction wind velocity V,^ 
is sufficiently large that saltation occurs for all of the particles in the size range from 
100 pm to 1000 pm. Various assessments of the threshold velocity or threshold friction 
velocity on Mars have been made (e.g., in the Viking Mars Engineering Model and in 
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refs. 6,7, and 8). In figure III -A -2 3 of the Viking Mars Engineering Model there is an 
indication that a wind of 40 m/s at a height of 1 m would be below the threshold value for 
all particles if Zq were 30 /im. In figure III-A-24 there is an indication that 40 m/s 
would be below the threshold velocity for particles larger than 600 pm in the mean model 
or 900 pm in the Max pg model atmosphere, even if Zg were as large as 1 mm. How- 
ever, these threshold values were calculated for "mean -surface -level" conditions (i.e., 
the mean equatorial radius of the planet). If threshold values are recalculated by the 
method used in the Mars Engineering Model for the mean atmosphere 4.5 km below the 
mean equatorial surface (representative of the proposed Viking 75 A prime and backup site 
and the Viking B prime site) and a value of 300 pm for Zg, it is found that 40 m/s at 1 m 
equals or exceeds the threshold velocity for particles of a 1000-pm or less diameter. 

The threshold -velocity values herein may be open to question. All of the values 
were based on data published in 1941 by Bagnold (ref. 1). More recent data by Chepil 
(ref. 9) in 1958 indicate lower values; the most recent data by Iversen, Greeley, Pollack, 
and White (ref. 10) in 1973 indicate higher values than those of Bagnold. Also, the effects 
of aerodynamic lift and the effects of cohesion of surface material were not considered in 
the derivation of any of these values. The most likely effect of these latter considerations 
would be a reduction of the threshold velocities for Mars. Thus, the assumption which is 
made here that all particles up to 1000 pm in diameter are involved in the saltation pro- 
cess is rather uncertain, but appears to be reasonable. Certainly it would be unconser- 
vative to assume that the larger particles were not involved. 

Computations 

Once the necessary input parameters are specified, the determination of apex 
height Zmax as a function of Wg for the assumed conditions on either Mars or Earth 
is a straightforward iterative process. Similarly, once particle velocities and scale 
height are also specified, the computation of mass flux kinetic energy, particle flux, the 
number of impacts for each particle class, the "effective" number of impacts, and sensor 
life is a straightforward process. However, since the exponential form of equation (3) is 
only an approximation of the actual variation of mass flux with height, the value of mass 
flux at a particular height depends on the way in which the fitting of the exponential is 
performed. In particular, the mass -flux value depends on the choice of Zj^ax either 
Mars or Earth used in applying equations (7) and (8). Larger values of Zj^^^x for either 
planet produce larger values of hj^/hg and larger values of mass flux on Mars in the 
region between 0.5 m and 2 m. The procedure adopted in the present report is: 

1. For each particle size and for each height z for which fluxes are to be computed 
for Mars, assume Zj^^x jyi equal to z. 
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2. From trajectory simulations find Wq corresponding to jyj. 

3. From equation (7) find Wq g corresponding to Wq 

4. From trajectory simulations find E corresponding to the Wq g from 

step 3. 

5. From equation (6) find hg. 

6. From equation (8) and h^, Zjjjax,M> ^max,E from the previous steps, 
find hj^. 

7. From equation (6) find Hjyj. (In calculations for sensitivity testing, the first 
seven steps are performed for only a single Mars height.) 

8. Using equations (2) and (3) find the value of mass flux Q for each height. This 
value is the equivalent mass flux for sand of a uniform size, and must be multiplied by the 
proportion of sand in each size class to find the mass flux of that size class for a sand 
with a wide range of particle size. From this point, the remaining flux calculations are 
straightforward. 


Sensitivity to Input Parameters 

The sensitivity of velocities computed in the trajectory simulations was tested by 
independently varying several of the input parameters over a range of values while holding 
the other parameter values constant at typical values. The ranges of parameters tested 
are: 

Zq = 30 jum to 300 jum 
k' = 3 mm to 1 cm 

p = 1.22 X 10"^ kg/m^ to 1.8 x 10”^ kg/m3 
jU = 1.6 X 10'^ to 10”'^ kg/m2-s 
j = 0.45 m/s to 2.5 m/s 

Except for the density change none of these parameters produced changes in the calcu- 
lated particle velocity greater than 1 m/s, or greater than about 3 percent of the velocity 
at heights above 0.5 m. The 50-percent change in density produced particle-velocity 
changes of up to 2 m/s, or up to 18 percent. 

The Zq value of 30 fxm corresponds to a smooth surface of sand with a predomi- 
nant size of 300 pm; Zq of 300 pm corresponds to measurements (not during saltation) 
of actual Earth desert surfaces. Measurements in wind tunnels during saltation are near 
the lower value. Values of Zq > 300 pm may occur on Mars if material other than the 
modeled sand is present. However, Chepil (ref. 11) and Iversen, Greeley, Pollack, and 
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White (ref. 10) show that the threshold velocity increases rapidly with increasing values 
of Zq] thus, it is unlikely that the larger sand particles would saltate. 

For Earth conditions Bagnold found k' to be about 3 mm. Zingg found k' to vary 
with particle diameter and to be equal to lOd (for an assumed mean diameter of Martian 
sand of 300 jum Zingg' s relationship also gives k' = 3 mm). In view of the lower gravity 
(3.7 m/s^) and possible greater saltation height of Mars, an additional value k' = 1 cm 
was used as a conservative value. 

The lower value of atmospheric density is the density at the mean equatorial surface , 
The higher density corresponds to an elevation 4.5 km below the mean equatorial surface 
at the prime and backup sites for the Viking 75 A mission, and the prime site for the 
B mission. Other possible landing sites are at higher elevations 

The viscosity values represent an arbitrary ±20-percent range about the estimated 
value of 1.8 x 10"^ kg/m-s^. 

Although the variation of particle velocity is small over the ranges of input param- 
eters considered herein, variation of mass flux may be substantial. Equation (2) indicates 
that height -integrated mass flux varies directly with p and as the cube of ; however, 
there are also indirect effects resulting from changes in H in equation (3) so that the 
total effects are best shown by a complete calculation as described herein. 

Combined Variation 

In order to illustrate the range of variation of particle velocity and mass flux, com- 
putations were made for both a "smooth" case and for a "rough" case. Parameters used 
for the two Mars cases, and for the single Earth case (used with both Mars cases) are 
given in table V. A 1-m wind speed of 40 m/s for Earth would result in values outside 
the range of values used in the formation of the empirical equations of Bagnold and Zingg. 
The Earth case is chosen to give values well within the ranges of experimental values of 
both Bagnold and Zingg. In the Mars cases parameters were chosen from the ranges 
previously used to give, in the smooth case, the highest wind speeds below 1 m and, in 
the rough case, to give the largest value of V/. In both cases a fixed 1-m wind speed 
of 40 m/s was employed. Velocities for the rough case are given in table III. Ratios of 
particle velocities for the two cases are given in table VI. The largest differences in 
particle velocities at heights where the kinetic energy is sufficient for sensor breakage 
(0.8 m and above) are about 5 percent compared with about 3 percent for variation of indi- 
vidual parameters. In table VII velocity ratios for the high-density and low-density rough 
cases are shown. Differences of up to 18 percent appear. 

To assess the effects of scale heights, mass -flux computations were performed for 
two different ratios of Mars scale height to Earth scale height for each case. For each 
case the total mass flux (for all particle sizes) was computed with the ratio of scale heights 
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for each particle size held constant with height, rather than varying with height as de- 
scribed in the "Computations" section. The two ratios were those computed for Mars 
heights of 0.1 and 1.6 m. The four resulting plots of mass flux plotted against height 
are shown in figure 9. Not only are the differences between the smooth and rough cases 
very large, but even larger differences are found between values computed using scale - 
height ratios determined for Mars heights of 0.1 and 1.6 m. 

In figure 10 computed mass -flux profiles are shown for the same wind-profile param- 
eters for density values of 0.12 kg/m^ and 0.18 kg/m^. Both of these mass -flux profiles 
represent the rough case described above and scale -height ratios computed using a Mars 
height of 1.6 m. As might be expected the calculated mass flux is smaller in the lower - 
density atmosphere. Although the variation here is smaller than that produced by changes 
in the wind -profile parameters, there are still order -of -magnitude changes in the mass 
flux at the larger heights. 

The particle -velocity calculations varied only a few percent over the range of Mars 
parameters considered. In contrast, the niass-flux model is quite sensitive to changes in 
the Mars parameters assumed. 


RESULTS AND DISCUSSION 

Because of the sensitivity of the mass -flux model to the variation of assumed Mars 
parameters, it would be difficult and perhaps meaningless to select a single average or 
typical set of Mars parameters for sensor expected-life calculations. Instead, attention 
is focused on the worst reasonably likely conditions, namely, the rough case with the 
higher -density value described herein. However, rather than using the scale -height 
ratios computed for 1.6 m (which produced the larger mass flux in the sensitivity tests), 
separate scale -height ratio values are computed for each Mars height as well as for each 
particle -size class, since this procedure gives the flux profile with the best fit at that 
particular height. 

Various flux values calculated by the methods described herein are presented in 
tables Vin, IX, and X for heights of 0.8, 1.2, and 1.6 m, respectively. In these tables P 
IS the percentage of particles of a given size class in the assumed surface material and 
equivalent mass flux is the total mass flux which would exist if all the particles were of 
the given size class; the mass flux of that size class is, then, the product of the equiva- 
lent mass flux and the percent P of particles of that size class. The mass flux is also 
equal to the product of the percent mass flux and the total mass flux. Similar relations 
hold for kinetic -energy flux, particle flux, and concentration. The number of hits per day 
for each particle size includes all the particles; the total number of hits per day and of 
effective hits per day is calculated only for particles which exceed the critical value of 
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15 pij. ("Hits per day" refers to the number of hits in a total period, not necessarily 
continuous, of 24 hours of saltation.) 

At 0.8 m, particles smaller than 700 jum do not attain the critical energy value of 
15 /zj. At 1.2 and 1.6 m, particles smaller than 400 fim do not reach 15 [ij. Thus, as 
suggested earlier, possible deficiencies in modeling the scale heights of particles smaller 
than about 125 jim do not affect the validity of the model for the sensor -breakage study. 

The dominance of the large particles at the greater heights is clearly apparent in the per- 
centage figures. 

The highest value of kinetic energy at 0.8 m (17 jitj) is only slightly in excess of the 
critical value of 15 iiJ , below which no breakage occurred; this value is also within the 
range where test subjects survived thousands of impacts. The highest value at 1.6 m 
(43 nJ) is triple the minimum value for breakage. This value coincides with 43 pj of the 
1200-/am particles in the test series. 

It IS also of interest, for example, to examine the characteristics of the complete 
flux profiles resulting from the combination of the profiles for individual particle sizes 
in order to determine whether this combination profile agrees with the assumption of an 
exponential variation of mass flux with height. Profiles of mass flux, concentration, and 
kinetic -energy flux are shown in figures 11, 12, and 13. The linear combination of indi- 
vidual exponential mass -flux profiles with scale heights which are not constant with height 
could not be expected to produce an exactly exponential variation with height. However, 
the validity of the model used depends on the assumption that an exponential function is 
a good approximation to the variation of mass flux with height. The mass -flux variation in 
figure 11 is a reasonable approximation to an exponential (which would appear as a straight 
line on this semilog plot), especially at the greater heights, although not as good an approx- 
imation as the curves in figures 9 and 10 where the scale height for each particle size was 
held constant over the height range. 

The concentration of sand (the mass of sand per unit volume of space) shown in 
figure 12 and the kinetic -energy flux shown in figure 13 both exhibit roughly exponential 
variation and show the dominance of the mass -flux variation over the particle -velocity 
variation. 

Table XI summarizes the application of the flux model to the problem of breakage 
of the hot -film wind sensors of the Meteorology Instrumentation System of the Viking 75 
Mars Lander. These data apply to the rough, high-density case. Expected sensor life- 
times for the smooth and low-density cases (not shown) are measured in thousands of 
years. The most critical level is 0.8 m, the lowest height at which the critical value 
of 15 /iJ occurs. Here the model yields a rate of 0.8 particle impacts per 24 hours of 
saltation with a 1-m wind of 40 m/s. Based on the assumption of a total of 7 days (from 
the Mars Engineering Model) of saltation per 90 mission days, and of a mimmum of 1000 
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effective impacts required for breakage, the expected lifetime of a sensor is 40 Earth 
years (=20 Martian years). This time period appears to be sufficient to assure a remote 
possibility of breakage, during either the 60-day primary mission or the proposed extended 
mission of one Martian year. The likelihood of breakage may be greater during an ex- 
tended mission than during the primary mission since landings are scheduled during the 
season when minimum winds are expected at the landing sites; however, the assumption 
of a total of 7 days of saltation during 90 mission days is considered to be sufficiently 
conservative to include year-round occurrences. 

At levels below 0.8 meter the number of impacts increases. However, the kinetic 
energy drops below the critical value of 15 pj and the present model does not predict 
breakage at these levels. 

At higher levels the number of impacts decreases rapidly while the kinetic energy 
increases. At the highest expected sensor level (1.6 m) the kinetic energy reaches the 
highest value encountered in the test data. Thus, values for greater heights are based on 
extrapolation of the test data but the continuation of the trend toward increasing expected 
lifetimes seems to be clear. 


CONCLUDING REMARKS 

A model for estimating the parameters of the flux of sand on Mars has been devel- 
oped using the wind-tunnel measurements of Bagnold and Zingg. The model is not sup- 
ported either by rigorous theoretical derivation or by experimental observations under 
realistic Martian conditions, but the assumptions involved appear to be both reasonable 
and conservative. Questions are raised concerning the validity of the model for particles 
which may be in partial suspension (<125 jam). These questions do not affect use of the 
model for studies of breakage of this wind sensor where such particles are shown to be 
ineffective. However, use of the model is restricted in material -erosion studies where 
such particles have been shown to be important. 

Application of this model to the sensor -breakage problem indicates that the expected 
lifetime of the wind sensors is adequate both for the primary Viking 75 mission and for a 
proposed extended mission with a duration of one Martian year. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
September 10, 1975 
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TABLE I. - FRACTURE CONDITIONS OF TEST SAMPLES 


Test 

sample 

Particle 
diameter, jim 

Particle 
velocity, m/s 

Kinetic 
energy, pj 

Number 

tested 

Number of 
collisions 

Range 

Mean 

Uncoated 

991 to 1397 

1200 

6 

43 

4 

25 to 50 

glass 

417 to 495 

450 

14 

13 

4 

250 to 350 

rod 

208 to 246 

225 

39 

12 

2 

2000 to 3000 

Prototype 

246 to 351 

300 

28 

15 

1 

20 000 


208 to 246 

225 

49 

19 

1 

24 000 

Flight 

type 

417 to 495 

450 

17 

18 

1 

^2100 


^The first 1000 collisions occurred at the subcritical speed of 12 m/s (9 /uJ), 


TABLE II.- SIZE DISTRIBUTION 

Percentage by mass of particles with mass nearest 
to the mass of a particle of given diameter 
^ (corresponds to dune -sand distribution of the 
Viking Mars Engineering Model) 


Diameter, 

jum 

Percentage 

100 

8 

200 

17 

300 

44 

400 

13 

500 

5 

600 

4 

700 

3 

800 

3 

900 

2 

1000 

1 
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TABLE in. - PARTICLE VELOCITY, m/s 
[computed using parameters for high-density rough case (table 
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4. 1 

3 . 5 

0.50 

32. 4 

•1 . o 

17.2 

13.0 
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7. 7 

C t< 
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■*1 cr 

. J 
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13. 1 
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11.6 

9.7 
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7 . Cl 

6 i 

1. 10 

j*' y ■ y 
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2 3 . 9 
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15.0 
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c- 

1 J 

b 5 

1.20 
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19.6 
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12. 8 
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9. 1 

7 . 

b. 9 
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16.3 

13.3 
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9. 5 

c ■ 3 

7 

1 . 40 
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3 . 1 

yy ■ w 
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9 . 9 

. 6 
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TABLE IV.- PARTICLE KINETIC ENERGY, juJ 
{computed using parameters for high-density rough case (table V)] 


100 

200 

300 

400 

500 

600 

700 

300 

900 

1000 

2-M 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

0 . 1 0 

0 

1 
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2 

d. 

Zf 

2 

2 


0. 20 

0 

~J 

c. 

Zf 

"3 

4 

4 

4 

4 

4 

4 

0. 30 

1 

"J 

4 

5 

S 

6 

6 

6 

6 

6 

0 . 40 

1 

5 

5 

6 

1 

8 

8 

8 

9 

9 

0 , 50 

1 

"5, 

6 

8 

9 

10 

10 

11 

1 i 

11 

0 . 60 

1 

4 

1 

9 

li 

12 

13 

1 3 

14 

14 

6, 70 

1 

4 

7 

10 

12 

14 

15 

16 

16. 

17 

0 . SO 

1 

4 

8 

12 

14 

16 

17 

18 

19 

20 

0. 90 

1 

ej 

9 

13 

16 

18 

20 

21 



1 , 00 

1 

S 

10 

14 

18 

20 

22 

24 

25 

36 

1, 10 

1 

5 

11 

16 

20 

22 

25 


28 

39 

1 , 20 

1 

b 

11 

17 

21 

25 

27 

39 

31 

“1 

1 . 30 

i 

b 

12 

18 

2 1 

27 

30 

2' 3 

34 

76 

1 . 40 

1 

s 

13 

19 

25 

29 

32 

35 

'• 7 

39 

1,50 

1 

ti 

14 

21 

26 

51 

35 


ki 

43 

1 , 60 

1 

1 

14 

22 

28 

33 

37 

41 

44 

46 

1 . 70 

1 

i' 

15 

U -• 

30 

35 

40 

43 

47 

49 

1.30 

1 

1 

15 

24 

31 

37 

42 

46 

50 

C ~i 

J S' 

1 . 90 

1 

1 

16 

25 

I 5 Z\ 

39 

45 

4'^ 

53 

56 

2 , 00 
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TABLE V.- PARAMETERS USED IN CALCULATIONS 



Mars high -density 
rough case 

Mars low -density 
rough case 

Mars high -density 
smooth case 

Earth 

p 

0.018 kg/m3 

0.0^22 kg/m^ 

0.018 kg/m^ 

1.225 kg/m^ 

p 

1.1 X 10-5 kg/m-s 

1.1 X 10"5 kg/m-s 

1.1 X 10"5 kg/m-s 

1.7894 X 10-5 kg/m-s 

g 

3.7 m/s2 

3.7 m/s2 

3.7 m/s2 

9.80665 m/s 

^0 

300 (im 

300 jim 

30 fxm 

30 jLtm 

k’ 

1 cm 

1 cm 

3 mm 

3 mm 

V at 1 m ... 

40 m/s 

40 m/s 

40 m/s 

11 25 m/s 

V*.t 

1.24 m/s 

1.24 m/s 

1.24 m/s 

0.19 m/s 

Vt 

10.87 m/s 

10.87 m/s 

14.28 m/s 

2.187 m/s 

V ' 



2.53 m/s 

2.53 m/s 

1.771 m/s 

0.624 m/s 


TABLE VI.- VELOCITY RATIO 


Ratio of velocity for high -density rough case to 
_ velocity for High -density smooth case (table V)_ 


100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

z,m 

0.89 

0.80 

0.73 

0.73 

0.74 

0.76 

0.74 

0.74 

0.75 

0.75 

0.10 

0.86 

0.83 

0.84 

0 83 

0.82 

0.83 

0.83 

0.82 

0.83 

0.82 

0.20 

0.91 

0.90 

0.87 

0.87 

0.87 

0.86 

0.87 

0.87 

0.87 

0.86 

0.30 

0.93 

0.91 

0.90 

0.90 

0.90 

0 89 

0.89 

0.89 

0.89 

0.89 

0.40 

0.95 

0.93 

0.92 

0.92 

0.91 

0.91 

0.91 

0.91 

0.91 

0.91 

0.50 

0.96 

0.95 

0.94 

0.93 

0.93 

0.93 

0.93 

0.93 

0.93 

0.92 

0.60 

0.97 

0.96 

0 95 

0.94 

0.94 

0.94 

0.94 

0.94 

0.94 

0.94 

0.70 

0.98 

0.97 

0.96 

0.96 

0.95 

0 95 

0.95 

0.95 

0.95 

0.95 

0.80 

0.99 

0 97 

0.97 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

0.90 

0.99 

0.98 

0.98 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

1.00 

1.00 

0.99 

0.98 

0.98 

0.98 

0.97 

0.98 

0.97 

0.97 

0.97 

1.10 

1.00 

0.99 

0.99 

0.98 

0.98 

0.98 

0.98 

0.98 

0.98 

0 98 

1.20 

1.01 

1 00 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.98 

0.98 

1.30 

1.01 

1.00 

1.00 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

1.40 

1.01 

1 00 

1.00 

1.00 

1.00 

0.99 

1.00 

1.00 

0 99 

0.99 

1.50 

1.02 

1.01 

1.01 

1 00 

1.00 

1.00 

1 00 

1.00 

1.00 

1.00 

1.60 

1.02 

1.01 

1.01 

1.01 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.70 

1.02 

1.01 

1.01 

1.01 

1 01 

1.01 

1.01 

1.01 

1 01 

1.01 

1.80 

1.02 

1.02 

1.02 

1.01 

1.01 

1.01 

1.01 

1.01 

1.01 

1.01 

1.90 

1.03 

1.02 

1.02 

1.01 

1.01 

1.01 

1.01 

1.01 

1.01 

1.01 

2.00 
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TABLE Vn.- VELOCITY RATIO 




Ratio of velocity for high-density rough case to 






velocity for low -density rough case 

(table V) 




D= 100 
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300 
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700 
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TABLE Vm.- FLUX VALUES AT 0.8 METER 
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TABLE IX.- FLUX VALUES AT 1.2 METERS 
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TABLE X.- FLUX VALUES AT 1.6 METERS 
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TABLE XI.- SUMMARY OF FLUX VALUES AND SENSOR LIFETIMES 
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Particle diameter , d , urn 

Figure 2.- Results of sensor -breakage tests. Symbols indicate particle velocity required 
to break the test subjects. Circles represent uncoated glass rods; squares represent 
prototype sensors; and the triangle, a flight -type sensor. Diagonal lines are constant 
values of kinetic energy. 
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Figure 9.- Mass flux as a function of height for the high -density rough and smooth cases described 
in table V for scale heights based on Mars height of 0.1 m and 1.6 m. 
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